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Abstract 

We study the excitonic effects on the second-order nonlinear optical properties of semi-spherical quantum dots 
considering, on the same footing, the confinement potential of the electron-hole pair and the Coulomb interaction 
between them. The exciton is confined in a semi-spherical geometry by means of a three-dimensional semi- 
parabolic potential. We calculate the optical rectification and second harmonic generation coefficients for two 
different values of the confinement frequency based on the numerically computed energies and wavefunctions of 
the exciton. We present the results as a function of the incident photon energy for GaAs/AIGaAs quantum dots 
ranging from few nanometers to tens of nanometers. We find that the second-order nonlinear coefficients exhibit 
not only a blue-shift of the order of meV but also a change of intensity compared with the results obtained 
ignoring the Coulomb interaction in the so-called strong-confinement limit. 



Introduction 

Nonlinear optical properties of semiconductor quantum 
dots have attracted considerable interest due to their 
several potential applications [1-4]. In particular, sec- 
ond-order nonlinear optical properties, such as non- 
linear optical rectification (OR) and second harmonic 
generation (SHG), have received special theoretical [5-8] 
and experimental [9,10] attention due to their magni- 
tudes being stronger than those of high-order ones, 
making them the first nonlinear optical effects experi- 
mentally observed. 

The confinement of carriers provided by a quantum 
dot is well described by a parabolic potential when only 
the lowest excited states of the carriers are considered. 
However, self-assembled quantum dots growth in the 
laboratory usually exhibit asymmetrical shapes that 
ensure the generation of nonlinear optical effects. In 
order to model these asymmetries, an asymmetrical 
potential is required. 

Recently, several authors [5,6,8] studied the effects of 
an exciton on the second-order nonlinear properties in 
one-dimensional semi-parabolic quantum dots. Using 
analytical approximate results, they showed that the 
excitonic effects enhance significantly the OR and SHG 
coefficients. They used the so-called strong-confinement 
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limit, ignoring in this way the Coulomb interaction 
between electron and hole because of the quantum dot 
dimensions are smaller than the effective Bohr radius, 
and finding that the excitonic effect reduces itself to an 
effective-mass model in one particle scheme. 

In this study we find eigenenergies and eigenstates of 
an exciton in a semi-spherical quantum dot solving the 
corresponding three-dimensional Schrodinger equation 
using a finite elements method and taking into account 
both the confinement and Coulomb potentials of the 
electron-hole pair. We present the OR and SHG coeffi- 
cients as a function of the incident photon energy with 
and without Coulomb potential. Our results show that 
energy and intensity of the peaks in the second-order 
nonlinear optical coefficients change when Coulomb 
interaction is introduced. 

This article is organized as follows. In "Theory" sec- 
tion, we present the characteristic quantities of the har- 
monic and Coulomb potentials, and the definitions of 
the weak- and strong-confinement limits in terms of 
these parameters. In addition, we present the analytical 
expressions for the optical nonlinearities, such as OR 
and SHG, obtained by the density matrix formalism. In 
"Results" section, we show the OR and SHG coefficients 
with and without Coulomb interaction as a function of 
the incident photon energy for two quantum dot sizes. 
We also give account of the changes presented by the 
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second-order nonlinear coefficients focusing in the role 
played by the Coulomb interaction. Conclusions are 
summarized in final section. 

Theory 

The effective-mass Hamiltonian for the electron-hole 
pair in the three-dimensional quantum dot reads [11] 



H 



_Ee_ 

2m* 



+ V(r e ) + 



2m* 



+ V(r h ) 



47T£|r e 



r h | 



(1) 



where m* and m£ are the effective masses of the elec- 
tron and the hole, respectively, s is the background 
dielectric constant and V (r,) is the three-dimensional 
semi-parabolic potential that we define as 
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The angle 6 is the usual polar angle in spherical coor- 
dinates, and co 0 the oscillator frequency considered in 
this study the same for the electron and the hole. The 
potential defined in Equation 2 confines the exciton in 
the upper half of a sphere, i.e., the quantum dot has a 
semi-spherical shape. 

Hamiltonian (1) can be separated in terms of center- 
of-mass and relative coordinates, respectively, 



H = 
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where M = m* + ra£ is the total mass, and 
[A = m^m^/M is the reduced mass. The center-of-mass 
and relative position coordinates are defined as usual, 



R: 



m*r e + m*r h 
M 



r h , 



(4) 



with the corresponding momenta P 
ihV T in terms of p e and p h , 



ih V R and p 
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The explicit separability of the center-of-mass and rela- 
tive coordinates in Equation 3 lead to the following total 
envelope wave function and total energy for the system: 



^(r c ,rh) = 4>(R)0(r), 



E = Er + E T . 



(6) 



(7) 



The center-of-mass part of Hamiltonian (3) is a three- 
dimensional semi-parabolic oscillator that can be solved 
analytically. Therefore, the problem is now reduced to 
solve the relative motion Hamiltonian: 



2fi 



47ts\r\ 



(8) 



Hamiltonian (8) has been solved analytically in two 
limiting cases (strong and weak confinement) for one- 
dimensional quantum dots. The eigenfunctions and 
eigenvalues are presented in references [5] and [8]. In 
one-dimensional case, the confinement potential also 
imposes constraints to spatial coordinates, resulting in a 
hydrogen-like (asymmetric-harmonic) reduced particle 
Hamiltonian for weak (strong) limit. 

The harmonic potential in Equation 8 defines both the 
size L of the quantum dot, 



h 

flCO 0 ' 



(9) 



and the energy quanta hco 0 due to confinement, which 
is related to L by 



h(L>o = -. 



fiL 2 



(10) 



On the other hand, the Coulomb potential defines the 
effective Bohr radius a$ and the effective Rydberg energy 
7£* of the electron-hole interaction, 



Ansh 
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The strong-confinement limit is established when 
L Uq, or equivalently hco 0 ^> 7^*, and the weak- 
confinement limit when L ^> aj, or hcoo <^ 7Z*. 

The second-order nonlinear optical coefficients can be 
obtained by density matrix approach and perturbation 
expansion method [12,13]. The expression for the OR 
coefficient, within a two-level system approach, is given 
by [5,6] 



(2) A e cr s 2 
Xo = 4—72 ^oi^oi 
£ 0 n 



l^oi 



T 2 



(ODQl + CO) + 



where e is the electron charge, a s is the density of 
electrons in the quantum dot, T 1 is the longitudinal 
relaxation time, T 2 is the transverse relaxation time, and 



Moi = I ifo\r\fi) \ 



(14) 
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For the resonance condition 

hco ^ Ei — Eq, 
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there is a peak intensity given by {II T x « II T 2 « 
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The SHG coefficient in a three level system is [8] 

e 3 N 
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where N is the density of carriers in the quantum dot, 
Eg = Ei - E J} r m = r 2 o = r 0 are the relaxtion rates, and 



Mij = \ [fi\r\fj)\, 



(20) 



Under the double resonance condition, i.e., hco « E 10 
E 20 /2, the intensity of the peak is given by 



(2) 2 3 N/zoi/xi2/x 2 o 
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and its energy by 
£20 + 2Eio 



(21) 



(22) 



Results 

In this study, the results are presented for a GaAs/ 
AlGaAs structure. We have used the following para- 
meters in the calculations: m* = 0.067m 0 , ^ = 0.09m 0 
(m 0 is the mass of a free electron) [4], T x = 1 ps, T 2 = 
0.2 ps [12], g s = 5 x 10 24 m 3 [5], £ = 12.53, T 0 = 1/ 
0.14ps ~\N=3 x 10 16 cm" 3 [8]. 

In Figure 1, we plot the characteristic lengths and 
energies for the confined particle in a GaAs/AlGaAs 
quantum dot as a function of the confinement frequency 
co 0 . cIq and 7£* are independent on co 0 because they are 
related to the Coulomb potential. In Figure la, we can 
see that the lengths L and are of the same order of 
magnitude for a confinement frequency around co 0 = 1 x 
10 13 s 1 . In Figure lb, we observe that also hco 0 and 
7^* show similar values around co 0 = 1 x 10 13 s" 1 . For 
this reason, we conclude that, in this frequency range, 
neither the strong-confinement limit nor the weak limit 





Figure 1 Characteristic (a) lengths and (b) energies for the 
confined particle in a GaAs/AIGaAs quantum dot as a function 
of the confinement frequency. The red (black) lines correspond to 
L and f)co 0 {a^ and 7£*), respectively. 



can be assumed because both interactions, harmonic 
and Coulomb, are important. Therefore, we propose a 
numerical technique to calculate eigenenergies and 
eigenstates of Hamiltonian (8), considering the harmo- 
nic and Coulomb potentials. 

With the aim of exploring the nonlinear behavior at 
higher frequencies, i.e., when the quantum dot size is 
smaller than in the previous case, we choose co 0 = 2 x 
10 14 s" 1 , in which the quantity L is less than cIq, or hco 0 
is greater than TZ* t differing in both cases by one order 
of magnitude as can be seen in Figure 1. Because of this 
difference, several authors [5,6,8] used the strong-con- 
finement limit as a satisfactory approximation in the 
case of small quantum dots. Accordingly with Equation 
9, the frequencies co 0 = 1 x 10 13 s" 1 and co 0 = 2 x 10 14 
s" 1 define a quantum dot size of L = 17 A nm and L = 
3.9 nm, respectively. This means that our results are sui- 
table for the current quantum dot sizes that range from 
few nanometers to tens of nanometers. 

We find numerically eigenenergies and eigenstates of 
Hamiltonian (8) by a finite elements method for the two 
frequencies mentioned above. We have used the soft- 
ware COM-SOL Multiphysics, which offers the possibi- 
lity of defining a geometry, in this case the upper half of 
a sphere, and to solve the time-independent Schrodinger 
equation with appropriate boundary conditions. 
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The terms involving quantum states and energies in 
Equations 13 and 19 are found using the eigenstates and 
eigenenergies previously calculated. The OR and SHG 
coefficients are shown in Figures 2 and 3, respectively. 
Figures 2a and 3a correspond to co 0 = 1 x 10 13 s" 1 , and 
Figures 2b and 3b to co 0 = 2 x 10 14 s" 1 . In each figure, 
we present the corresponding nonlinear optical coeffi- 
cient considering excitonic effects with and without 
Coulomb interaction. For comparative purposes, we also 
present the case without excitonic effects, i.e., when 
only one electron exists in the quantum dot. 

Figures 2 and 3 reproduce the reported results quite 
well [5,8] regarding the enhancement of the nonlinear 
optical coefficients due to the quantum confinement of 
the exciton. However, the same figures show addition- 
ally a significant blue-shift of the OR and SHG peaks 
when both the harmonic and Coulomb potentials are 
taken into account. In Tables 1 and 2, we present the 
eigenenergies of the exciton and peak energies of the 
OR and SHG coefficients with and without Coulomb 
interaction for the two frequencies under study. The 
peak energies are estimated by Equations 17 and 22 for 
the OR and SHG coefficients, respectively. 




2JQ 340 250 260 270 2S0 200 300 
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Figure 2 The OR coefficient as a functions of the incident 
photon energy t)co for (a) co 0 = 1 x 10 13 s" 1 and (b) co 0 = 2 x 
10 14 s"\ considering excitonic effects with (red line) and 
without Coulomb (black line)interaction. The blue line 
corresponds to the case without excitonic effects. 

v / 
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Figure 3 The SHG coefficient as a functions of the incident 
photon energy hco for (a) co 0 = 1 x 10 13 s" 1 and (b) co 0 = 2 x 
10 14 s"\ considering excitonic effects with (red line) and 
without Coulomb (black line)interaction. The blue line 
corresponds to the case without excitonic effects. 



We can see from Tables 1 and 2 that the eigenener- 
gies obtained with Coulomb inter-action are smaller 
than those obtained without that interaction. The expla- 
nation to this fact is that there is an attractive Coulomb 
potential between the electron-hole pair that implies a 
reduction of the eigenenergies for the exciton. However, 
the eigenenergies are affected in different ways depend- 
ing on the quantum state. For example, for the ground 
state co 0 = 1 x 10 13 s~\ Table 1, we have an energy dif- 
ference of 5.237 meV, while for the first and second 
excited states the differences are of 3.472 and 2.784 



Table 1 Eigenenergies of the exciton and peak energies 
of the OR and SHG coefficients with and without 
Coulomb interaction for co 0 = 1 x 10 13 s" 1 



Energy (meV) 


With Coulomb 


Without Coulomb 


Diff. 




11.218 


16.455 


5.237 




26.147 


29.619 


3.472 


E 2 


40.000 


42.784 


2.784 


OR peak energy 


14.929 


13.164 


1.765 


SHG peak energy 


14.660 


13.164 


1.496 
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Table 2 Eigenenergies of the exciton and peak energies 
of the OR and SHG coefficients with and without 
Coulomb interaction for co 0 = 2 x 10 14 s" 1 



Energy (meV) 


With Coulomb 


Without Coulomb 


Diff. 




306.59 


329.10 


22.51 


E ] 


577.06 


592.39 


15.33 


E 2 


843.31 


855.67 


12.36 


OR peak energy 


270.47 


263.28 


7.19 


SHG peak energy 


269.41 


263.28 


6.13 



meV, respectively. We have a similar situation for co 0 - 
2 x 10 14 s" 1 , Table 2. This is because the mean spatial 
separation between the electron and the hole increases, 
and therefore the Coulomb interaction decreases, as the 
energy of the quantum state increases. The final result 
is a blue-shift of the OR and SHG peaks of the order of 
meV for both co 0 = 1 x 10 13 s" 1 and co 0 = 2 x 10 14 s" 1 . 

In addition, the OR and SHG coefficients exhibit dif- 
ferent peak intensities depending on the consideration 
of the Coulomb interaction, as it can be seen in Figures 
2 and 3. This fact originates from the modification of 
the dipole matrix elements defined in Equationsl4, 15, 
and 20 when the Coulomb interaction is considered. 
According to Equation 18, the peak intensity of OR 
coefficient depends essentially on the product /Xq^oi, 
while for SHG coefficient, Equation 21, the peak inten- 
sity depends on fi 01 [i 12 ^ 2 o Tables 3 and 4 show the 
values of these dipole matrix element products with and 
without Coulomb interaction for the two frequencies 
considered. As one can see, the product /Xq^oi is greater 
with Coulomb interaction than without it for both con- 
finement frequencies. Therefore, in Figure 2a, b, the OR 
intensity is higher in the former case than in the later 
one. In the case of SHG coefficient, the product 
P20 is smaller with Coulomb than without that interac- 
tion. This fact makes the SHG intensity smaller in the 
former case as can be seen in Figure 3a, b. 

Conclusions 

Contrary to the assumption that Coulomb interaction 
can be neglected when the quantum dot dimensions are 
smaller than the effective Bohr radius, we show that this 
interaction affects the excitonic effects of the second- 
order nonlinear optical properties of semi-spherical 
quantum dots. We find that Coulomb interaction 

Table 3 Dipole matrix element products of the OR and 
SHG coefficients with and without Coulomb interaction 



for co 0 = 1 


x 10 13 s" 1 






Coefficient 


[nm 3 ] 


With Coulomb 


Without Coulomb 


OR 




1365 


1222 


SHG 


M 0] M U M 20 


1237 


1635 



Table 4 Dipole matrix element products of the OR and 
SHG coefficients with and without Coulomb interaction 



for co 0 = 2 x 10 14 s 1 


Coefficient 


[nm 3 ] 


With Coulomb 


Without Coulomb 


OR 




16.95 


14.03 


SHG 


^011^1 2^20 


17.65 


18.55 



manifests itself in a blue-shift of the energy peaks of the 
order of several meV in the studied spectra. These results 
were found for two quantum dot sizes, in the first one 
the characteristic quantities of the harmonic and Cou- 
lomb potentials are equals, and in the second one they 
differ by one order of magnitude. This means that the 
Coulomb interaction plays an important role even when 
the quantum dot sizes are smaller than the effective Bohr 
radius. 

Therefore, we encourage experimentalists to carry out 
measurements of second-order optical nonlinearities in 
asymmetrical quantum dots with the aim of to detect 
the magnitude of this effect. 
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OR: optical rectification; SHG: second harmonic generation. 
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